Multi-scale clustering in granular surface flows 



D. Bonamy, F. Daviaud, L. Laurent, M. Bonetti, and J. P. Bouchaud 

Service de Physique de I'Etat Condense, CEA Saclay, 91191 Gif sur Yvette, France 

We investigate steady granular surface flows in a rotating drum and demonstrate the existence 
of rigid clusters of grains embedded in the flowing layer. These clusters are fractal and their size 
is power-law distributed from the grain size scale up to the thickness of the flowing layer. The 
implications of the absence of a characteristic length scale on available theoretical models of dense 
granular flows are discussed. Finally, we suggest a possible explanation of the difference between 
velocity profiles observed in surface flows and in flows down a rough inclined plane. 
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Granular materials share properties with both usual 
liquids and solids. They can form an inclined free sur- 
face without flowing, but when the angle of the free sur- 
face exceeds some threshold value, an avalanche occurs. 
Global behaviours can be described by models derived 
from fluid mechanics jj], || or non linear physics j|| Q. 
However, some experimental results remain far from be- 
ing understood. Detailed measurements of the mean den- 
sity profile and the mean tangential velocity profile ob- 
tained in three-dimensional flow (3D) by NMR || and 
in two-dimensional flows (2D) by direct image analy- 
sis |], [jj , show strong evidence that the relation between 
stress and strain is non local: (i) the velocity gradient is 
found to be constant in the flowing layer whereas momen- 
tum balance predicts a linear variation of the shear stress 
with depth || ; (ii) the velocity gradient presents a differ- 
ent scaling with the depth for dense granular flows down 
a rough inclined plane Jl(j, [ll], [l2| indicating the non local 
influence of boundary conditions on internal rheology in- 
side the flowing layer; (iii) the velocity gradient does not 
vanish at the free surface at variance with typical fluids. 
Very recent models have tried to account for some non 
local effects || ||, [|. We report here the first experi- 
mental evidence of rigid clusters of grains embedded in 
the flow and characterize their geometrical and statisti- 
cal properties. Although clustering instabilities driven by 
the inelasticity of grain collision are well known in granu- 
lar gases Jl3|, Q , they have never been observed in dense 
surface flows. We find that these clusters are fractal and 
their size is power-law distributed from the microscopic 
scale - the diameter of a grain - up to the macroscopic 
scale - the flowing layer thickness. Therefore, no charac- 
teristic correlation length can be defined in the flowing 
layer. Almost 50% of flowing beads belong to these ob- 
jects. We discuss recently proposed non local models in 
the light of these experimental observations, and propose 
a argument to explain the differences observed between 
granular surface flows and granular flow down a rough 
inclined plane. 

Experimental setup: The experimental setup is illus- 
trated in Fig. |l|a. It consists in a duralumin rotating 
drum of diameter Do = 45 cm and variable gap, half- 
filled with steel beads of diameter d = 3 ± 0.05 mm. The 
use of weakly dissipative beads of millimetric size pre- 



vent any clustering due to material-specific effects such 
as capillary or electrostatic effects. Two drum thickness 
are used: 

(A) a gap of 7 mm so that a quasi-2D packing is ob- 
tained but with a local 3D microscopic disordered 
structure. The fast camera allows then to track the 
beads actually seen through the transparent side 
wall of the tumbler: around 40% of the beads are 
hidden and cannot be tracked. 

(B) a gap of 3 mm so that a pure 2D packing is obtained 
and all the beads can be tracked. Additional metal- 
lic obstacles (bound to the inside of the drum) have 
been added in the static part of the packing to pre- 
vent any 2D ordering effects which would induce 
non generic effects. The main drawback is that the 
measurements can be only made during the time 
when the obstacles are in the bottom part of the 
drum. 

In the recorded region, located at the center of the drum, 
the flowing layer thickness R and the mean angle 9 of 
the flow can be changed by varying the rotation speed 
n. Regimes obtained for f2 varying from 1 rpm to 8 rpm 
are investigated. On this range of rotating speeds, sur- 
face flows are steady and inertial effect are negligible (the 
Froude number Fr = flDo/2g < 0.01 where g is the grav- 
ity constant). 

The beads are lighted via a continuous halogen lamp. 
Sequences of 200 frames are recorded via a fast camera 
at a sampling rate of 1 kHz. The recorded region size in 
pixel is 480 x 234, one pixel corresponding to 0.227 mm. 
Frame processing allows us to obtain the position of the 
center of mass of the beads seen through the transparent 
side wall of the tumbler. Since the image of a single 
bead is made up of about 23 ± 3 pixels (depending on the 
distance to the side wall) , the errors on the determination 
of bead location is about 150 /zm [jl^]. Tracking each bead 
on ten successive frames allows us to evaluate its velocity 
averaged on 10 ms (see Fig. |l|b). The mean angle 9 of 
the flow is calculated from the averaging of the velocity 
of all the beads on the whole sequence of 200 frames. The 
images are then divided into layers 1 bead diameter wide 
parallel to the flow. The mean tangential velocity V m (z) 
of the layer at depth z is then defined as the average of 



2 



0°00„ ° a Q Q Q^QCJ Q Q 



° ° o o°q ra Q _ Q 
° °o < 

m r To O O O o ° O O O „ O O o _ ra Q O S 

b o o o o o°s o o, G o o ° o ro^^.VJ A; 

3 O 8 - 




OQ^ °o e0 o° 0O 8„° 




(a) 



Fast camera 
(up to 1kHz) 



Lamp 





1500 


° £1=8 rpm 


o 


E 


1000 


□ £1=3 rpm 


o 



— <? 


E, 






s 

□ 


E 

> 


500 


*f2=1 rpm 











' (c) 



-80-60-40-20 20 40 
z (mm) 



FIG. 1: (a) Sketch of the experimental setup, (b) A typical 
instanteous velocity field obtained for quasi-2D flows (geom- 
etry A) with Q — 8 rpm in the rectangular box of (a), (c) 
Velocity profile measured at the center of the drum for 3 dif- 
ferent rotation speed: Q = 1 rpm, f2 = 3 rpm and = 8 rpm. 
The arrows indicates the end of each profile. The precision 
on each point is better than the point size. Velocity profiles 
are linear in the flowing layer with a gradient independent of 
Q, i.e. independent of both the local slope 6 and the flowing 
layer thickness R. 




FIG. 2: (a) Typical instantaneous field of the orientation of 
the beads velocity fluctuations obtained for a quasi-2D flow 
with Q — 8 rpm. The orginal instantaneous velocity field 
is represented in gray, (b) Distribution function P(8) of the 
angle S between the orientation of the velocity fluctuations 
of two beads in contact. The straight line corresponds to a 
uniform fonction. 



the velocity of all the beads of the sequence of 200 frames 
whose center of mass is inside the layer. 

For both geometries (A) and (B), the granular surface 
flow indeed presents a clear linear velocity profile (see 
Fig. [l]c for the quasi 2D case) with a velocity gradient 
7 = dV m /dz independent of both R and 0. In the follow- 
ing, we neglect beads with velocities less than 100 mm/s 
since these are assumed to belong to the static bed which 
solidly rotates with the drum. 

Velocity fluctuations in quasi-2D and pure 2D flows: 
In order to investigate possible collective effects, we 
study first the spatial correlations of the instantaneous 



velocity field for both geometry (A) and (B). Calling 
~V(x,z,t) the velocity of a bead of the frame t located 
at the (x,z) coordinate in (e x ,e z ), where the unit vec- 
tor e x (resp. e 2 ) is parallel (reps, perpendicular) to 
the mean flow, the bead fluctuation velocity is defined 
as: V(x, z, t) = \{x, z, t) — V m (z)e x . Since the orienta- 
tion of the fluctuation n = V/V is not correlated to the 
depth z (at variance with the fluctuation amplitude V), 
we focussed on the instantaneous field made up of these 
orientations. 

One of these instantaneous fields is represented on 
Fig.^a and reveals aggregates of strongly correlated 
beads. These correlations can be quantified by looking at 
the distribution of the angle S between the orientation 
and hj of two beads i and j in contact Q (Fig. ||b) . This 
distribution reveals that the velocity fluctuations of two 
beads in contact tend to have correlated orientations. To 
isolate the correlated aggregates, we have used the fol- 
lowing criteria: Two beads i and j belong to the same 
cluster whenever (1) they are in contact and (2) the angle 
between and is smaller than a given value arbitrary 
chosen equal to 60 °. This value corresponds to the point 
where the experimental orientation distribution crosses 
the value of the uniform distribution (see Fig. ||b). Val- 
ues ranging from 30° up to 90 ° have also been tried and 
do not modify the following conclusions. The Figs. ||a,b 
show typical frames of the individual clusters for geome- 
try (A) respectively for = 8 rpm and fi = 1 rpm. Note 
that the size of the largest ones are of order of the flow- 
ing layer thickness R. We also plot the number of beads 
N in a cluster against the cluster radius of gyration R g 
(Fig. |c). The power law scaling observed in both ge- 
ometries (A) and (B) is indicative of a fractal structure 
with an apparent fractal dimension equal respectively to 
d { f A) = 1.6 ±0.1 and df = 1.9 ±0.1 for quasi 2D and in 
pure 2D flows. This quantifies the observation that the 
clusters are more open and stringy in the quasi 2D geome- 
try. More interestingly, the distribution of the number of 
beads N in a cluster is a power-law N~ a in both geome- 
tries (A) and (B), with an exponent equal respectively to 
= 2.9 ±0.1 and = 2. ±0.1 for quasi-2D and 2D 
flows (see Fig. [p for the quasi-2D case). In both cases, 
the exponent is independent of Q, i.e. independent of 
both R and 9. This indicates that no typical correlation 
length scale can be defined to characterize the spatial 
correlations of the velocity field. 

Volume fraction fluctuations in pure 2D flows: In order 
to find out the physical origin of these correlated clusters, 
measurements on the local volume fraction fluctuations 
were performed on the 2D stacking (geometry (B)). A 
Voronoi tesselation can then be used to define the local 
volume fraction v associated with each bead. Let us note 
that this technique assumes that all beads position are 
known and consequently cannot be applied for the quasi- 
2D packing of geometry (A) . The relative fluctuations of 
v are very small, a few percent (see Fig. ^a), but sufficient 
to alter significantly the beads behaviour. According to 
Reynolds' dilatancy concept fL7|, an assembly of rigid 
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FIG. 3: Clusters of beads with correlated velocity fuctuation 
orientation in quasi 2D flows (geometry (A)), (a) and (b): 
Typical frames of the clusters respectively for £2 = 8 rpm and 
= 1 rpm (see text for details en the cluster definition). 
The position of the beads belonging to the flowing layer are 
represented with a black circle. The gray arrows represent 
the raw velocity field. Beads belonging to the same clus- 
ter are drawn in the same color (c): cluster size TV plotted 
against their radius of gyration R g . The errorbar indicated 
the standard deviation on R g of clusters with a given N. A 
power law behaviour (straight line) indicates fractal struc- 
tures of dimension d 



(A) 
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1.6 ± 0.1. (d): Probability density 
function of the size (in number of beads N) of these clus- 
ters. Data have been shifted for clarity. From bottom to top, 
£2 = 1 rpm, £7 = 2 rpm, Q — 3 rpm, Q — 4 rpm, £7 = 6 rpm 
and £7 = 8 rpm. For these different £7, i.e. for different 8 and 
R, the distribution decreases as a power-law (straight line) 
with the same exponent a^ A ' = 2.9 ±0.1. 



particles cannot deform whenever its volume fraction is 
higher than a given threshold v c . For pure 2D packing 
of monodisperse beads, v c — 7r/4. Clusters of beads in 
contact with v > v c have been isolated in the flowing 
layer (see fig. f|d,e,f). Almost 50% of beads of the flowing 
layer belong to one of these solid clusters. As was found 
for the clusters defined using velocity correlation, the size 
of largest ones is of order of R. The size distribution is 
again a power law with an exponent ah ~ 1.5, weakly 
dependent of £7 (see Fig. ||b). This exponent is slightly 
smaller than the one obtained from by analysis of velocity 
fluctuation. It should be noted that although there is 
some overlap between both types of clusters they are not 
identical. The fractal dimension of these 'dense' clusters 
is found to be d v = 1.7 ± 0.1. Their mean aspect ratio, 
calculated as the ratio of the two principal axis of their 
inertia tensor, does not depend on their size and is found 
to be close to 2. They preferentially orient in the flow 
direction (see Fig. |^c). Their life time, defined as the 
duration during which their size is larger than the half of 
their maximal size, is of order of \J d/g ~ 10 -2 s, where 
g is the gravity. 

Discussion: These solid clusters can be compared to 
the collective objects postulated in different non local 




FIG. 4: Presence of rigid clusters in the flowing layer of pure 
2D flows (geometry (B)). (a) Averaged volume fraction profile 
u m (z). (b) Probability density function of the size (in number 
of beads N) of the " rigid" clusters (see text for details) . Data 
have been shifted for clarity. From bottom to top, £7 = 2 rpm, 
£7 = 7 rpm, £7 — 10 rpm and £7 = 13 rpm. They are power 
law distributed with an exponent close to ai B ' ~ 1.5 weakly 
dependent of £7. (c) Polar diagram of the probability density 
of cluster direction, showing a preferential orientation parallel 
to the flow. Frames (d) (e) and (f ) show a typical sequence of 
these clusters. Beads belonging to the static phase appear in 
gray. Beads belonging to the same cluster are drawn in the 
same color. The clusters are tracked on the successive frames. 



models for granular flows. Mills et al. (D postulate the 
existence of ID transient solid chains in the flowing layer, 
well separated from each other. In their approach, that 
describes granular flow down a rough inclined plane, a 
grain colliding with one of these chains transmits its 
momentum throughout the whole chain. On the other 
hand, for surface flows, Rajchenbach || suggests that 
each shock impact is imparted to the whole substrate 
beneath and thus that the momentum is 'dissipated' into 
the whole packing. This description has been recently 
extended by Andreotti and Douady || who account for 
the possible trapping of flowing grains in the bumps of 
the static bed. The localization of the flow within a layer 
of finite thickness and the linear velocity profile are well 
reproduced in their approach. However their model does 
not account for the experimental dependence of the ve- 
locity gradient in the thickness R and the local angle 9. 
The collective objects that we observe are neither chain- 
like nor the whole packing, but rather widely distributed 
fractal clusters involving 50% of the beads of the flowing 
layer. No typical length scale can be defined. This last 
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point may be responsible for the selection of a velocity 
gradient independent of both R and 9, i.e. independent 
of the shear stress. 

The observation of the successive clusterised frames 
shows that the largest clusters are emitted by the static 
phase and die either by fragmenting into smaller ones or 
by sticking back to the static phase. Differences between 
granular flows down a rough inclined plane and gran- 
ular surface flows can then be rationalized as follows: 
clustering results from the competition between inelas- 
tic multiple collisions, that tends to aggregate grains to- 
gether jl3|, [uj, and shear, that erodes clusters. For flows 
down an inclined plane, these two effects lead to clusters 
of a typical size, whereas in surface flows, the static bed 
plays the role of a cluster reservoir: its erosion by the 
flowing grains can generate very large size clusters that 
then split and cascade into smaller and smaller ones (this 
mechanism might explain why we observe power-law dis- 
tributions). In this view, the cluster size distribution, 
and a fortiori the velocity gradient, depend crucially on 
the boundary conditions. 

The absence of a characteristic correlation length to 



describe the locally 'jammed' clusters of beads is quite in- 
teresting. First, this contrasts with an assumption made 
in most local and non local models [§ |, § §, where 
the transition between the 'solid' phase and the 'liquid' 
phase is supposed to occur over a well defined length 
scale. The existence of multi-scale rigid clusters indicates 
that the flowing phase is actually critical, and suggests 
the proximity of a continuous 'jamming' transition, of the 
type proposed in |Q (see also (l9) for a related discus- 
sion). In this respect, it is interesting to note that similar 
power-law distributed clusters of strongly correlated mo- 
tion have recently been observed in a colloid close to the 
glass transition pC| . 
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